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Abstract We present grids of stellar models and their as- 
sociated oscillation frequencies that have been used by the 
CoRoT Seismology Working Group during the scientific pre- 
paration of the CoRoT mission. The stellar models have been 
calculated with the CESAM stellar internal structure and evo- 
lution code while the oscillation frequencies have been ob- 
tained from the CESAM models by means of the ADIPLS adi- 
abatic oscillation programme. The grids cover a range of 
masses, chemical compositions and evolutionary stages cor- 
responding to those of the CoRoT primary targets. The stel- 
lar models and oscillation frequencies are available on line 
through the Evolution and Seismic Tools Activity (ESTA) 
web site. 
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E> 1 Introduction 

, During the scientific preparation of the CoRoT mission, a ref- 
erence frame in the H-R diagram had to be provided for the 
CoRoT potential target stars. In particular, the critical selec- 
tion of the primary targets and the census and characterisa- 
tion of the stellar content in the fields surrounding the po- 
tential primary targets have required to locate all considered 
stars (or stellar systems) in the H-R diagram and eventually 
to estimate their mass, evolutionary stage and expected os- 
cillation spectra. For that purpose we have calculated the in- 
ternal structure and evolution of stars in a range of masses, 
chemical compositions and evolution stages of interest for 
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the modelling of CoRoT targets (iMichel and Baglin'. '2005'), 
We have used the CE SAlvP] stell ar evolution code (MoreJ, 



1 19971 : iMorel and Lebretoiin2007h . We have selected several 
models along each evolution sequence and we have calcu- 
lated their oscillation frequencies by mean s of th e ADIPL^ 
programmes (IChristensen-Dalsgaardlll982i 120071) . 

In this paper we briefly present these grids of models 
and oscillation frequencies. Section |2] presents the CESAM 
numerical code and the input physics and initial parameters 
used to calculate the grids while Section |3] presents the os- 
cillation data obtained from the ADIPLS adiabatic oscillation 
code. The present models and associated oscillation frequen- 
cies have been made available on the ESTA web siteB 



2 CoRoT/CESAM stellar models 

2.1 Input physics and numerical aspects 

The CESAM stellar evolution code has been written and de- 
veloped by R Morel and collaborators since t he late 1980's. 

The ge neral descrip tion of t he cod e is given bv lMorel and LebretonI 
(l2007h but see also' Morell (Il997h . CESAM is a public code. 



The source is available on-line at the web site - http : //www . obs-nice . f r/c 
- together with a short guide of directions for use and a com- 
prehensive document giving a complete description of the 
numerics and of the physics implemented. All the models 
presented here were calculated in 2003 with the 4.2 version 
of CESAM written in Fortran77. We point out that the present 
version of CESAM, CESAM2k, was re -programmed in For- 
tran95; it has also been updated with new input physics and 
has been en riched with the inclus ion of new physical pro- 
cesses (Mo rel and LebretonI l2007l) . 

The models have be initialised on the zero age main se- 
quence (ZAMS) and the evolution has been followed up to 

' Code d' Evolution Stellaire Adaptatif et Modulaire 
^ Aarhus Adiabatic Pulsation Package), available at 
http : //astro . phy s . au . dk/ j cd/ adipack . n 
^ http://www.astro.up.pt/corot/models/ and 

also on the Paris Observatory web site at 

|http : //wwwTisr . obspm . f r/^lebreton/Modeles/CESAM_CDRDT . html] 
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the beginning of the red giant branch. The models have typi- 
cally 600 mesh points in the interior while the atmosphere is 
restored on a grid of typically 100 grid points. During hydro- 
gen burning the time step is adjusted according to the relative 
changes of the hydrogen abundances. In the present calcula- 
tions it takes from 165 to 420 time steps to evolve from the 
ZAMS to the ascent of the red giant branch, depending on the 
mass. 

All the models have been calculated with the same given 
set of standard input physics. Here, we have neglected mi- 
croscopic and turbulent diffusion processes as well as rota- 
tion or magnetic fields. The input physics are the following: 



(1 19931) corresponding to a solar metals to hydrogen ratio 
{Z/X)q = 0.0245 where Z and X are the abundances in 
mass fraction of heavy elements and hydrogen respec- 
tively. In the nuclear reaction network the initial abun- 
dance of each chemical species is split between its iso- 
topes according to the isotopic ratio of nuclides. The ini- 
tial amount of is assumed to be converted in ^He. 
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- Equation of State . We used the so-called CEF F analytical 

equation of state dChristensen-Dalsgaard an d Dappen, 1992b 

which is an improved version of the EFF dEggle ton. Faulkner, c'2.0 
Il973h equation of state including the Coulomb correc- 
tions to pressure in the so-called Debye-Huckel approx- 
imation. ^ 

- Opac ities. We used the OPAL95 opacity tables dlglesias and R( 
Il99d ) complemented at low temperatures {T ^ \(f'K) 
by the Alexander and Ferguson (1994) tables. The metal 
mixture of the opacity tables is fixed to the initial mix- 
ture of the model (see below). 

- Nuclear reaction rates. We used the basic pp and CNO 
reaction networks up to the 0{p,ay'^N reaction. In 
the present models the CESAM code takes ^Li, ^Be and 

at equilibrium. The nuc lear reaction rates are from 
ICaughla n and Fowled(ll988h . Weak screening is assumed 
under Salpeter's formulation (|1954|) . 

- Convection and oversh ooting. We use t he cla ssical mix- 
ing length treatment of Bohm-Vitense ('1958') under the 
formulation of Henyey, Vardya, and Bodenheimer ( 1965}) 
taking into account the optical thickness of the convec- 
tive bubble. To find the location of the onset of convec- 
tion, we used Schwarzschild's criterion. In models calcu- 
lated with overshooting, the convective core is extended 
on a distance /qv = CLo\ x min(//p,7?cc) where Oov is the 
overshooting parameter, the pressure scale height and 
/?cc the radius of the convective core. The core is mixed 
in the region corresponding to both the convective and 
overshooting regions and in the overshooting region the 
temperature gradient is taken to be equal to the adiabatic 
gradient. 

- Atmosphere. Eddington's grey T{t) law is used for the 
atmosphere calculation. The hydrostatic equation is in- 
tegrated in the atmosphere, starting at the optical depth 
T = 10~^ and the connection with the envelope is made 
at T = 1 to ensure the vali dity of the diffusion approx- 
imation dMorel et fl/.LIl994i) . At this level, we insure the 
continuity of the variables and of their first derivatives 
with space. The radius of the star is taken to be the bolo- 
metric radius, i.e. the radius at the level where the local 
temperature equals the effective temperature (t = 2/3 
for the Eddington's law). 

- Initial abundances of the elements and heavy elements 
mixture. All models were calculated with the classical 
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Fig. 1 Grid Gl: evolutionary sequences in the H-R diagram for mod- 
els with solar metallicity [Fe/H] = 0.0 coiTesponding to Z/X = 0.0245, 
helium abundance Y = 0.2672, no overshooting. 



2.2 Initial parameters of the models 

The calibration of a solar model in luminosity and radius 
with the input physics given above leads to an initial abun- 
dance of helium Yq = 0.2672, and to a solar mixing-length 
parameter ttiviLT,© = 163. The corresponding initial hydro- 
gen and metal abundances are respectively =0.7153 and 
Zq = 0.0175. 

We have calculated models for two values of the [Fe /H] 
ratio: [Fe/H] = 0.0 (solar) and —0.1 (metal deficient). The 
metallicity of the models Z has been obtained from the rela- 
tion [Fe/H] = log(Z/Z) - log(Z/Z)Q. The individual abun- 
dances of metals have been scaled on the solar ones. 

For models with non solar metallicity, we have derived 
the helium abundance in mass fraction Y according to the 
relation Y = Y^+Z{AY /AZ), where AY / AZ is the helium 



1.3andyp = 



to metal enrichment ratio and Y-p is the primordial helium 
abundance. We iodkAY /AZ ^ (Ffg 
0.245 (see Olive and Skillmanir2004 . 

All the models were calculated with a mixing-length pa- 
rameter equal to the solar one. Finally, for solar metallicity 

we ca lculated two grids of models, one without overshooting 

GN93 solar mixture of heavy elements from lGrevesse and I'tadab ne including overshooting with the value ttov = 0.15. 
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We have calculated models with masses in the range cov- 
ered by the CoRoT targets: 1 .0, 1 .2, 1 .4, 1 .6, 1 .8, 2.0, 2.2, 2.4, 
3.0,4.0, 5.0 and 10.0 Mg. 




Fig. 2 Grid G2: evolutionary sequences in the H-R diagram for mod- 
els with solar metallicity [Fe/H] = 0.0 corresponding to ZjX = 0.0245, 
helium abundance Y = Q.lifll and overshooting Oov = 0.15 



2.3 Model grids 

We present three grids of models. Grid Gl consists in evolu- 
tionary sequences for the 12 values of the stellar mass listed 
above (between 1.0 and 10.0 Mq), a solar metallicity and 
no overshooting. This grid covers the domain of mass of 
main CoRoT targets, including solar-like oscillators and 5 
Scuti and j3 Cephei pulsators. In addition to this standard 
grid, we have calculated two grids restricted to solar-type 
and 5 Scuti stars, i.e. for the 8 lower values of masses listed 
above, between 1.0 and 2.4 M©. Grid G2 consists in evo- 
lutionary sequences with solar metallicity and overshooting 
(aov = 0. 15) and grid G3 consists in evolutionary sequences 
with [Fe/H] = —0.1 and no overshooting. Figures [l] |2] and 
[3] display the evolutionary sequences in the H-R diagram 
for the 3 grids of CESAM models available on the ESTA web 
site0. For each evolutionary sequence, we also provide on 
the web site a file in which the properties of the models 
are displayed: mass, logarithm of the effective temperature 
log Teff, bolometric magnitude, logarithm of the luminosity 
log(L/L0), logarithm of the surface gravity log g and age. 



3 ADIPLS oscillation frequencies 

Using the ADIPLS stellar oscillation program, we have cal- 
culated the oscillation frequencies for a selection of models 
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Fig. 3 Grid G3: evolutionary sequences in the H-R diagram for mod- 
els with sub-solar metallicity [Fe/H] = —0.1 coiTesponding to Z/X = 
0.0194, helium abundance Y = 0.2632 and no overshooting 

of the G2 grid (i.e. the grid corresponding to models with 
solar metallicity [Fe/H] = 0.0, ZjX = 0.0245, helium abun- 
dance Y = 0.2672 and overshooting aov = 0.15). In relation 
with CoRoT targets, we have chosen two values of the mass 
in the G2 grid: M = 1.2 Mq which is representative of so- 
lar like oscillators and M = 1.8 Mq which is representa- 
tive of 8 Scuti pulsators. For the M = 1 .2 Mq sequence, we 
have selected models along the evolutionary track by steps 
of 100 Myr, except on the hook at the end of the main- 
sequence where steps are of 10 Myr. For the M = 1.8 Mq 
sequence, we have selected models by steps of 50 Myr on the 
main-sequence and by steps of 10 Myr from the red point to 
the bottom of the red giant branch. 

For each mass we provide two archives on the ESTA web 
site, a Oi and a Fi archive, where i=l and 2 for M = 1 .2 and 
M = 1.8 Mq respectively. 

The Oi archive contains ASCII files, with a .osc exten- 
sion, each corresponding to a selected model on the evolu- 
tionary sequence. The .osc files can be read by means of the 
programme tools available in the ADIPLS package and de- 
scribed on the ESTA web sit^ They list the global properties 
of the considered model (name, mass, age, initial composi- 
tion, luminosity, photospheric radius, etc.) as well as a wide 
range of model variables at each mesh point (distance r to 
the centre, mass inside the sphere of radius r, pressure, den- 
sity, temperature, chemical composition, opacity and sev- 
eral other physical parameters including quantities of in- 
terest for the computation of adiabatic oscillations like the 
Brunt-Vaisala frequency). 

The Fi archive contains ASCII files, with a freq exten- 
sion, each corresponding to a selected model on the evolu- 
tionary sequence. The .freq files are the result of the fre- 
quency calculation by the ADIPLS program. They list the 



* |http : / /www . astro . up . pt/corot/models7 



^ at |http: //www. astro. up.pt /corot/ntools/modconv 
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mass, radius central pressure and density of the model and 
the frequency of oscillation, for each mode identified by its 
degree £ and order n. We point out that these files are not 
corrected for problems of numbering of the modes that oc- 
cur for modes of ^ = 1 in the region of the fundamental. 



4 Conclusion 

On the ESTA web site we have provided grids of stellar mod- 
els and oscillation frequencies that have been used by the 
CoRoT Seismology Working Group during the scientific prepa- 
ration of the CoRoT mission. The stellar models are provided 
for a range of masses between 1.0 and 10.0 M© and chemi- 
cal compositions [Fe/H]=0.0 and —1 .0 dex representative of 
the main CoRoT targets (solar-type stars, 5 Scuti stars and j3 
Cephei stars). Oscillation frequencies are provided all along 
the main-sequence track for a solar-type oscillator {M =1.2 
Mq) and a 5 Scuti pulsator (M = 1.8 Ma). Up to now these 
reference grids have been used to locate CoRoT potential tar- 
gets in the H -R diagram in the proc ess of target selection 
(see Fig. 1 in iLebreton et al '., 2007 ^) and to s tudy some 5 
Scuti candidates for CoRoT ( Poretti etalX. I2005h . 

On the ESTA web site, other grids of stellar models and 
oscillation frequencies can be found. These stellar models 
have been calculated either with the CESAM or the CLES 
code and the oscillation frequencies have been calculated ei- 
ther with the ADIPLS, POSC or LOSC code. All this material 
(codes a nd grids) is described in this volume by Lebreton et al\ 
d2007al) : lMontalban et a/1 (l2007a):' M"arques et a/.l ( l2007t) . We 
point out that the thorough comparisons of the stellar evolu- 
tionary codes and oscillations codes performed under ESTA 
have shown very good agreement between models calcu- 
lated by the CESAM and CLES codes whic h gives us con 



fidence in the use of the present grids (see 



l2007bl : lMontalban et a/.ir2007bl : lMova et al. 



Lebreton et all 
20071) . 
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